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ABSTRACT 

We report a detection (3.5 x 10 37 ± 5.6 x 10 36 ergs -1 ) of the optical coronal emission 
line [Fe x]A6374 and upper limits of four other coronal lines using high resolution VIMOS 
spectra centred on NGC 4696, the brightest cluster galaxy in the Centaurus cluster. Emission 
from these lines is indicative of gas at temperatures between 1 and 5 million K so traces the 
interstellar gas in NGC 4696. The rate of cooling derived from the upper limits is consistent 
with the cooling rate from X-ray observations (~10 M yr -1 ) however we detect twice the 
luminosity expected for [Fe X] A6374 emission, at 1 million K, our lowest temperature probe. 
We suggest this emission is due to the gas being heated rather than cooling out of the intr- 
acluster medium. We detect no coronal lines from [Ca XV], which are expected from the 5 
million K gas seen near the centre in X-rays with Chandra. Calcium is however likely to be 
depleted from the gas phase onto dust grains in the central regions of NGC 4696. 
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1 INTRODUCTION 

The X-ray emitting gas at the centers of clusters of galaxies can 
have very short cooling times. Without a heating source this gas 
should cool and condense at rates of up to thousands of solar masses 
per year. However, high resolution X-ray spectroscopy has been 
unable to detect large quantities of gas cooling below one third of 
the cl uster virial temperature (for a review see IPeterson & Fabianl 
l2006h . Observed rates of star formation and quantities of cool 
and cold gas in brightest cluster galaxies (BCGs), the most mas- 
sive galaxies known, a re far lower than expected (see for example 
iMcNamara et al.l2006h . For the hot gas to be radiating but not cool- 
ing in the predicted quantities a heating mechanism is required. 

To determine how the heating and cooling balance in clus- 
ters of galaxies is maintained requires a multi- wavelength approach 
to probe the many phases of matter. The hot intra-cluster medium 
(ICM), seen in X-rays, is at a temperature of ~ 10 7 — 10 8 K. Cool, 
optical emission line gas, found to surround the majority of BCGs 
in cool core clusters is seen at temperatures of ~ 10 4 K. Cold 
molecular gas reservoirs have also been observed, however obser- 
vations of gas at intermediate temperatures (10 5 — 10 6 K) have thus 
far proved elusive. 

High-ionisation, collisionally-excited optical coronal lines are 
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emitted at temperatures between 10 5 — 10 6 K making them im- 
portant tracers of intermediate temperature gas. Optical emission 
lines have the advantage of being observable from the ground. The 
high spectral and spatial resolution in principle allow us to deter- 
mine and trace the velocity structure, and therefore examine the 
transport processes, in the gas. Emission lines from these plas- 
mas can also prov ide a direct means of measuring the mass flow 
rate JCowidfl98lh and therefore th e rate at which the ICM is con- 
densing. iGraney & Sarazml dl990h . ISarazin & Granevl dl99lh and 
IVoit et al.l 4l994T) have modelled coronal line emission in conditions 
appropriate to those of the cores of galaxy clusters. 

There have been several attempts to detect coronal line emis- 
sion from cooling gas in clusters of galaxie s, howeve r repor ts of 
a significant detection remain inconclusive. iHu et al ] Jl985h ob- 
served the core regions of 14 clusters of galaxies but were unable 
to detect the presence of the [Fe XIV 1 A5303 in any of their long 
slit spectra. Similarly iHeckman et al.l dl989h attempted detection 
of another coronal line, [Fe X]A6374, in a sample of nine suspected 
'cooling flow' clusters. Likewise they found no detectable emis- 
sion. 

The first rep orted detection of coronal emission was by 
lAnton et all dl99lh in Abell 1795. The authors claimed to have 
found a 13a detection of the [Fe X]A6374 line in the core of the 
cluster with an aperture centred o n the central 2" of the cD g alaxy. 
A later study of the same object bv lDonahue & Stockel \ 19941) failed 
to confirm this detection. The sample used by iDonahue & Stockel 
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Figure 1. Top left: X-ray temperature map of the centre of the Centaurus cluster IFabian et all (2005). The colour bar is in units of keV. Box 1 and 2 are 
situated at the regions of lowest X-ray temperature. Top right: X-ray metallicity map based on Fe abundance. The colour bar is in units of solar metalicity. The 
0.4 — 0.8 keV X-ray temperature contours are overlaid. Bottom: [N IIJA6583 emission in NGC 4696 overlaid with X-ray contours of 0.4 — 0.8 keV gas. All 
images are on the same scale (10" = 2.1 kpc). The black X indicates the centre of the radio emission. Box 1 is at position RA 12^48 m 49 s .4, Dec. —41° 18' 
48.9" (J2000) and box 2 is situated at RA 12^48 m 48 s .7, Dec. -41° 18' 38.6" (J2000). Both boxes are 5 by 5 fibres large corresponding to 3.4 arcseconds 
(0.7 kpc) across. 



contained five massive cool core clusters and in this same 
paper they report three more non-detections (A2199, 2A0335+096 
and A2597) and one marginal (3a) detection in PKS 0745-191. 
This detection has a lso failed to be confirmed in a later study by 
lYan& Cohen! (Il995h . 

The emission line nebulosities around NGC 1275, the BCG i n 
the Perseus cluster were studied bv lShields & Filippenkd Jl992n . 
They searched for [Fe X]A6374 emission in off-nuclear regions, 
away from the central AGN, but also obtained no detection. 

Up until now the tw o most recen t searc hes for coronal 
line emission have been bylYan&Cohen1 ( ll995h and lSabra et alj 
bOOCh . lYan & Cohen] dl995h searched for emission from [Fe 
X]A6374 in Abell 1795 and PKS 0745-191, the only two clusters 
where detections of coronal line emission had previously been re- 
ported. Contrary to these results they did not find any significant 
detection in either cluster but were able to place upper limits on the 
surface brightness a factor of 10 deeper. 

In a larger program st udying the properti es of the emission line 
nebulosity in NGC 1275. ISabra etail bOOCh searched for coronal 
lines over a wide wavelength range using Keck. Six coronal line 
transitions were covered in this range from [Ni XII]A4232 to [Ni 
XII]A6702. No coronal line emission was seen in their data. 

This work reports on deep, high spectral resolution, inte- 
gral field spectroscopy (IFS) observations of NGC 4696 (RA 
12^48 m 49 s , Dec. -41° 18' 40" J2000), the BCG in the Centaurus 
cluster, Abell 3526. The Centaurus cluster at redshift z = 0.0104 
is the second nearest example of a cool core cluster. The heating 



and cooling in this cluster is apparently very well balanced despite 
the short central cooling time of only 200 Myr. NGC 4696 houses 
a radio source, and multiple bubbles. These are accompanied by 
soft X-ray filaments and a sharp rise in the metal abundance in the 
cen tral 30 kpc, among the highest seen in any cluster, ~ twice so- 
lar dSanders & Fabianll200i IFabian et al.ll2005l : ISanders & Fabianl 
l2006ah . Centaurus also has the broadest range of X-ray tempera- 
tures seen, cont aining gas from 0.35 to 3.7 keV, over a factor of 10 
m temperature dSanders et alJ2008h . 

Crawford et al ] b005h presented images showing the exten- 
sive, filamentary Ha nebulosity surrounding NGC 4696. These 
share the morp hology of the soft X -ray filaments, and of a promi- 
nent dust lane dSparks et al.lll989D . The low X-ray temperatures, 
short cooling times and exceptionally high metallicities in this 
nearby cluster make it the ideal target for a deep search for coronal 
line emission. 

The observations are briefly described in § 2, method, analysis 
and limits on coronal line emission is given in § 3 and in § 4 we 
discuss the implications of these limits. We summarise our results 
and main conclusions in § 5. At the redshift of the Centaurus cluster 
(z = 0.0104, 44.3 Mpc) one arcsec corresponds to 0.210 kpc 
(throughout this paper we adopt Ho = 71 km s _1 Mpc -1 , Qm = 
0.27 and Qa — 0.73). Ab undances given in th i s pape r are relative 
to the solar metallicities of lAnders & Grevessel dl989h . 
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Figure 2. The best-fit BC03 SSP model (red) to the observed spectrum, 
from box 1, black between 5700 — 6700 A and the fit residuals. Regions 
where we expect object emission lines and regions of poor sky subtraction 
were masked out in the fit and are shown above in grey. 
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Figure 3. Top: The range of STARLIGHT fits over the [Fe X]A6374 
region, in box 1, determined by perturbing and re-fitting the observed 
spectrum. Bottom: The distribution of integrated flux values for 200 
STARLIGHT fits over the [Fe X] A6374 region. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations were made on 2009 March 27-30 using the Visible 
MultiObje ct Spectrograph (VI MO S) on the VLT in La Paranal, 
Chile (see lLeFevre et alJl2003l and IZanichelli et al.ll2005l for a de- 
scription of the VIMOS IFU and a discussion of data reduction 
techniques). We obtained High Resolution Orange (HRO) data us- 
ing VIMOS Integral Field Unit (IFU). We used the larger 0.67" 
fibres giving a field of view of 27" x 27" with the HR grism. 

We acquired 10.5 hours exposure centred on the inner region 
of NGC 4696 (12 h 48'49.3", -41 o 18'40"), the deepest observations 
thus far to try to detect emission from cooling hot gas in these ob- 
jects. We also took six 15 minute exposures using the same set up 
to image the entire galaxy including many sky fibres. 
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Figure 5. Spectrum from box 1, after subtraction of a smooth continuum, 
showing the fit to the [O IJA6300 object emission line and the scaled tem- 
plate for the [O IJA6363 emission (red). The blue Gaussian is the result of 
a fit to the remaining excess flux after subtraction of the oxygen line. The 
dashed spectrum is the original continuum subtracted data and the solid line 
shows the effect of the [O IJA6363 line subtraction. Here also two neigh- 
bouring sky lines are fit and subtracted. The continuum is allowed to vary 
from zero and is fit locally around the emission lines of interest. 

The data were reduced by the VIPG0 pipeline (see 
IScodeggio etaD 120051 for a description of VIPGI). The 3D- 
datacubes were combined and analysed with a set of IDL routines 
(R. Sharp, private communication). 

Instrument flexures in VIMOS are very dependent on rotator 
position and suffer from hysteresis (lAmico et aDH008). For this 
reason we took calibration frames after each observation block. 

Cosmic-ray rejection, final fibre-to-fibre transmission correc- 
tions, sky subtraction, correction for Galactic extinction and shift- 
ing to the object rest frame were performed outside of VIPGI using 
IDL routines. The data reduction procedure will be explained thor- 
oughly in a forthcoming paper, Canning et al. (in prep), detailing 
observations of the emission line nebulae surrounding NGC 4696. 



3 METHOD, ANALYSIS, DETECTIONS AND LIMITS 

Our HR Orange spectra cover the wavelengths of a number of 
coronal lines, specifically [Fe XIV]A5303, [Ca XV]A5445, [Ca 
XVJA5694, [Fe X]A6374 and [Ni XV]A6702. 

The expected location of the [Fe X]A6374 feature coincides 
with the red wing of the [O IJA6363 object emission so requires 
special attention (see section [3^2l ). None of our other expected coro- 
nal lines coincide with object emission lines or sky lines. 

3.1 Continuum subtraction 

When determining flux from possibly broad and faint emission 
lines, errors in continuum subtraction are likely to dominate. For 
this reason continuum subtraction was attempted in a number of 
ways. The first method we employ is to use a small sample of 
fibers from within the galaxy, in a region where no significant coro- 
nal emission is seen. A second me thod is to subtract a smoo thed 
continuum from the spectra (as in iDonahue & Stockel 1 19941) and 

1 VIPGI- VIMOS Interactive Pipeline Graphical Interface, obtained from 
http ://cosmos . iasf -milano . inaf . it/pandora/ 
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Figure 4. Top: Spectra from box 1 (black) and box 2 (red) without continuum subtraction. All sky lines except [O IJA6300 and [O IJA6363 sky emission have 
been subtracted in order to more clearly see object emission line features. Bottom: A blow-up of the expected [Fe XJA6374 region. Emission from cool gas in 
box 2 has a broader velocity width than box 1 in general except over the [O IJA6363 and [Fe XJA6374 region. This demonstrates the excess flux on the red 
wing of oxygen emission in box 1. This can not be accounted for by subtraction of a template from the object [O IJA6300 emission line. 



finally we model t he continuum with simpl e stellar populations 
(SSP) models from lBruzual & Charlotl J2003h . hearafter refered to 
as BC03. These methods are described below. 

1. In an initial search we binned the spectra based on X-ray 
temperature, the lowest temperature regions (box 1 and 2, Fig. [T} 
overlap with the temperature range at which some coronal lines are 
emitted, specifically the lines of [Ca XV] (5 million K). In box 1 we 
found evidence for a flux excess on the red wing of the [O I]A6363 
emission line, blueshifted slightly from the 6374 A where we would 
expect [Fe XJA6374. We found no significant excess in box 2. We 
then use the spectrum of box 2 to correct the continuum of box 1 
and search for significant flux excess in regions of expected coronal 
emission. 

2. A smooth continuum is fitted to the sky subtracted spectra 
with all regions where emission lines are expected and regions of 
poor sky subtraction masked out. The smooth continuum is then 
subtracted from those spectra. 

3 . The continuum shape is fit usi ng the STARLIGHT soft- 
ware dCid Fernandes et al]|2005l [2QQ9h and basefiles made up of 
240 BC03 models spanning 6 metallicities, covering the range 
0.005 — 2.5 Zq, and 40 ages, covering the range — 20 Gyr. Using 
this software requires that we re-grid our spectra to integer wave- 
lengths. STARLIGHT uses a Markov Chain approach to fitting the 
data. We determine the parameters to use by fitting the data many 
times and looking at the distribution of fits over the region of in- 
terest. The fit with the best x 2 value is taken as the continuum and 
subtracted from the spectrum. An example of a fit with emission 
lines masked out is shown in Fig. [2 STARLIGHT does not allow 
us to examine the uncertainty in the parameters from a fit and so we 
use a Monte Carlo approach to determine the error. We quantify the 
additional uncertainty introduced due to the continuum subtraction 
by perturbing our data points by a Gaussian random number with 
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Figure 6. The sky spectra extracted from fibres at the edge of the field ~ 35" 
from the centre of NGC 4696. The [O I] A6363 sky emission line is labelled 
as are the regions where we would expect the object [O I] emission lines and 
the [Fe X] emission. The spectrum has been shifted to the rest- wavelength 
of NGC 4696. 



mean and variance given by the data and associated poission error. 
These new perturbed spectra are then fit with the STARLIGHT 
software using the same parameters as before. We repeat this 200 
times for our box 1 spectrum, the range of fits over the [Fe X] A6374 
region can be seen in the top panel of Fig. [3] The lower panel of Fig. 
[3] shows the distribution of the integrated flux (counts-model) over 
the wavelength region 6355 — 6380 A. The Gaussian shape of this 
distribution implies we have repeated the fits enough times to prop- 
erly sample the range of models which could provide a good fit to 
the data. The 200 model continua are then re-gridded to the same 
wavelength grid as our original data. We estimate the additional 
uncertainty in each spaxel as the one sigma deviation in the range 
of models. This uncertainty is dependent on the error spectrum of 
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Figure 7. Top: Box 1 spectra after continuum subtraction using the con- 
tinuum from box 2 spectra. Middle: Box 1 spectra with smooth continuum 
subtracted. Bottom: Box 1 spectra with best fit STARLIGHT continuum 
model subtracted. In each plot the top panel shows the data over the wave- 
length region containing the [O IJA6363 object emission and probable [Fe 
X]A6374 object emission. The red line shows the scaled [O IJA6300 emis- 
sion used as a template to subtract the contribution of the [O IJA6363 line. 
The bottom panel, in each case, shows the residuals from a subtraction of 
the [O I] A6363 template. In each plot an excess of flux can be seen around a 
wavelength of 6368 A. The two emission lines that have been greyed out, to 
the left of the [O IJA6363 line, are sky lines. These were fit and subtracted 
from the spectrum before a fit to the remaining excess flux was performed. 
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Figure 8. Box 2 spectrum with a smooth continuum subtracted. There is no 
significant excess emission seen on the red wing of the object [O IJA6363 
emission line. 
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Figure 9. The STARLIGHT fit in the bluest region of our spectra. After 
subtraction of the model stellar continuum we are left with variations in the 
residual spectrum on the scale of 10 — 20 A. This is comparable to the width 
of feature we find on the red wing of the [O IJA6363 object emission line. 
For comparison the line shows the ±lcr width of a 475 kms -1 FWHM 
velocity width line at the rest wavelength of the [Fe XIVJA5303 coronal 
emission line. 



the entire wavelength range of the spectra being fit and as such will 
only be weakly correlated with the possionian error on each spaxel. 
We therefore estimate the total uncertainty as the quadrature ad- 
dition of the poisonion error and the uncertainty in the continuum 
model. 

An advantage of using synthetic spectra for the continuum 
subtraction is that we can ensure that there is no coronal line emis- 
sion in the fit, so we are not subtracting any signal. The stel- 
lar spectra are much more complicated than a smoothed profile 
and we would like to test whether any apparent excess flux found 
where we would expect coronal line emission could be due to a 
'bump' in the stellar spectra. However the spectral models assume 
the relative abundances are the same which may not be appropri- 
ate to NGC 4696. The stellar continuum will also vary across the 
field and in binning spectra over a large region of the galaxy we 
would expect the SSP model fit to deteriorate. For these regions fits 
were run seperately on all spectra i.e. we did not use the best fit 
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STARLIGHT model of box 1 to correct for the stellar continuum 
of any spectra other than box 1 . 

3.2 [Fe X]A6374 

The [Fe XJA6374 coronal emission feature coincides with the red 
wing of the [O I]A6363 emission and as such warrants special treat- 
ment (see Fig. [5J. Unfortunately the object [O I]A6300 emission 
also coincides with the sky [O IJA6363 emission feature. 

We simultaneously fit Gaussians to the lines of [O I], the [N II] 
doublet, the [S II] doublet, Ha and the [Fe X] emission lines. This 
is in order to allow for the errors in scaling and removing the [O 
I]A6363 object emission line. A single Gaussian for each emission 
line in the obtical nebulosity gives a good fit in the regions of box 
1 and 2 (see upper panel of Fig.HJ). 

We tie the redshift and velocity width of the object optical emi- 
sion line nebulae ([N II] doublet, [S II] doublet and Ha and [O I]) 
emission lines. The redshift and velocity width of the sky [O I] lines 
are also tied. The integrated flux of the [N II] doublet and that of 
the two [O I] lines can be tied to each other, the scaling in each cas e 
being dictated by atomic parameters JOsterbrock & Ferlandll2006h . 
The integrated flux of the [Fe X]A6374 emission line is fixed to be 
a positive value but the redshift, and velocity width are allowed to 
be free in the fit as there is no evidence yet for a spatial correspon- 
dence between the optical nebulosity and any coronal line emitting 
plasma. The continuum level for each line is fit in a region local to 
the emission line. An example of the fitting process with the fit to 
the object [O I] emission lines and the [Fe X]A6374 emission line 
is shown in Fig. [5] 



3.3 Detections and upper limits 

After continuum subtraction, detections and upper limits are estab- 
lishe d by fitting a Gaussian line pro file to the spectra using MP- 
FIT dMorelll978l : lMarkwarddl2009l) . The 90 per cent flux limits 
are found by increasing the integrated flux in the fit until A% 2 has 
grown by 2.7, these results are presented in columns 2-4 of table 
Q] The last three columns show, for two cooling rates and metal- 
licities, the predicted luminosity of the coronal l ines assuming they 
are th e product of cooling from the hot ICM (see lSarazin & Granevl 
Il99ll) . During the fit we fix the Gaussian normalisation to be posi- 
tive. When determining the upper limits we also constrain the red- 
shift and velocity width to be the same as that of the [Fe X]A6374 
emission line. The coronal lines probe a broad range of tempera- 
tures so it is perhaps not obvious that they should have an identical 
velocity. The continuum is fit in a region local to the emission lines. 

We report a 63a detection of excess flux at the expected wave- 
length of [Fe X]A6374 emission in a region 15" away from the 
nucleus (box 1, Fig. [TJ. The quoted significance of this feature is 
from a fit to the SSP model continuum subtracted spectrum. This 
method of continuum subtraction gave the least significant result. 
We do not see any evidence of this emission in a region of the 
galaxy where there is similar low temperature X-ray gas nearer to 
the core (box 2). Fig.|4]shows a comparison of the spectra in these 
two regions. The velocity width of the cool (10 4 K) gas near the 
core is larger than that of the gas farther away in all cases except 
the region around the [O I]A6363 object emission. Here the box 1 
spectra shows an excess of flux towards the red wing of the oxygen 
line, this is not seen in the box 2 spectra or in the line profile of the 
[O I]A6300 line in box 1. 

Fig. [6] shows the sky spectrum. This spectrum is formed from 



regions at the edge of the field of view where we have no detectable 
object emission. The regions where we would expect object [O 
I] A6300, [O I] A6363 and [Fe X]A6374 emission is indicated. There 
is no obvious sky emission in the region around the [O I] A6363 and 
[Fe X]A6374 lines and as the excess flux is seen as a very broad 
feature and it is not seen everywhere in the spectra it is unlikely to 
be confusion due to sky emission. 

To fit the excess flux requires a redshift different from that of 
the optical emission line nebulae, the line in box 1 is blueshifted 
with respect to the optical nebulosity. The line of sight velocity dif- 
ference between the redshift of the Ha nebulae and the best fits to 
the excess flux is ~ 230 km s" 1 . This is comparable to the spread 
of the line of sight velocities in the optical nebulosity itself. The re- 
sults from the three methods of continuum subtraction agree within 
the error (see Fig. [7}. In box 1 the integrated flux with subtraction 
of the continuum using a region where no significant excess flux 
is observed, after subtraction of the [O I]A6363 line, gives an inte- 
grated flux of 1.8 x 10" 16 ± 2 x 10" 17 erg cm" 2 s"\ subtraction 
of a smoothed continuum gives an integrated flux in the line of 
1.6 x 10" 16 ±2x 10" 17 erg cm 2 s 1 , in both these cases the er- 
ror is determined from the possionian error in the spectrum which is 
similar to the one sigma deviation after continuum subtraction of a 
nearby emission free region of the spectrum, and finally SSP model 
fitting gives 1.5 x 10" 16 ± 2.4 x 10" 17 erg cm s , including 
both the possionian error and the additional error in the continuum 
fit. Fig.[8]shows, for comparison, the same region with the template 
[O I]A6363 line subtracted. There is no significant excess emission 
seen to the red-wing of the [O I]A6363 line. 

In regions containing strong stellar continuum features where 
the background spectrum is highly structured (A < 6000 A) the er- 
ror after continuum subtraction is larger than the simple poissonian 
error. 

The noise features in these regions have a width similar to the 
emission feature we are looking for (10 — 20 A, see Fig. [9}. We 
estimate the noise in these regions in two ways. First we calculate 
the one sigma deviation from zero in an emission free region of 
the spectrum (after continuum subtraction) on the scale of the vari- 
ations in the continuum. Second we fit a Gaussian of fixed width 
equal to the variations in the continuum. We then step the Gaus- 
sian over the emission line free region pixel by pixel. Due to the 
fluctuations in the noise the Gaussian fit will sometimes have a 
positive normalisation and sometimes a negative one. The distri- 
bution of the Gaussian area values in this region provide a second 
mechanism for estimating the one sigma uncertainty. The errors 
determined in both fashions were consistent with each other. Fig. 
[lOl shows the upper limits overplotted on the continuum subtracted 
spectrum for the [Fe XIV]A5303, [Ca XV]A5445, [Ca XV]A5694, 
and [Ni XV]A6702 lines. 

3.4 Spatial distribution of the hot gas 

The detection of [Fe X]A6374 emission in box 1 but not in box 2 
implies there may be some 'dumpiness' to the spatial distribution 
of the intermediate temperature gas. A small change in temperature 
equates to a vast change in emissivity (see Ferland et al. in prep), 
so this clumpy appearance may be indicative of small variations in 
temperature of the gas across the field of view. 

Using IFS observations we can trace the morphology of the 
regions where there appears to be an excess of flux towards the 
red wing of the [O I]A6363 line. Fig. Qj] shows the spatial distri- 
bution of this excess flux across t he field of view, o verlaid with 
the contours of radio emission (see lTavlor et al We binned 
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Figure 10. The 90 per cent upper limits overplotted on the continuum subtracted spectra from box 1 for the lines of [Fe XIVJA5303, [Ca XVJA5445, [Ca 
XV] A5694 and [Ni XVJA6702. Here the continum was subtracted using SSP models fit with the STARLIGHT software package. The redshift and velocity 
width of the emission lines are constrained to be the same as that of the detected [Fe XJA6374 line. 



Line 


Box l. a (10 37 ) 


Luminosity ( erg s 
Box 2 a (10 37 ) 


20 arcsec 2 . b (10 37 ) 


Z = Z Q (10 37 ) 


Model ( erg s 1 ) 
Z = 0.5Z© (10 37 ) 


Z = Z Q (10 37 ) 


[Fexiv]A5303 


<4.5 


<4.0 


<11.5 


3.1 


3.1 


10.3 


[Caxv]A5445 


<1.0 


<0.9 


<9.1 


1.4 


1.2 


4.7 


[Caxv]A5694 


<3.8 


<1.4 


<3.2 


2.1 


1.8 


7.0 


[Fex]A6374 


3.5±0.56 


<7.0 


25.8±8.3 


4.1 


4.2 


13.7 


[Ni xv]A6702 


<1.1 


<0.2 


<1.7 


0.2 


0.2 


0.6 



Table 1. The 90 per cent (2.7a) upper limits on the luminosity and detections of [Fe XJA6374 emission with one sigma errors in the regions specified. The 
luminosity quoted as detections used the continuum subtraction technique which gave the lowest integrated flux, in all cases this was when using the SSP 
model fits. Where the results are upper l imits and not clear detect ions the continuum subtraction method yielding the largest integrated flux was taken. The 
model approximate luminosities are from lSarazin & Granevl Jl99lh and scale simply with the cooling rate. Here they are given for a cooling rate of 3 (Columns 
5 and 6) and 10 Mq yr _1 (Column 7). These limits are determined by spectral fitting of the X-ray spectrum ( Sanders et al. 2008). The sizes of the regions 
in square arc seconds (kpc) are a 3.4 arcsec 2 (0.7 kpc 2 ) and b 20.1 arcsec 2 (4.2 kpc 2 ). These limits have been corrected for galactic extinction but not for 
the intrinsic reddening of NGC 4696. The upper limits on [Ca XV] A5694 imply this is deficient by a factor of 2 in the inner 20 arcseonds compared with the 
luminosity predicted from the X-ray cooling rate within this region. 



the spectra in 5 by 5 regions, this size binning was used as a com- 
promise between spatial resolution, signal to noise and goodness 
of fit to the stellar continuum. We then fit each region with SSP 
models as decribed in section ITT1 The [O I]A6363 object emission 
line was then fit and removed and the remaining 'counts' between 
6363 — 6380 A were summed (no emission line shape was assumed 
here) to create the map of relative flux, seen in the bottom panel of 
Fig. [TT] As an illustration we show, in Fig. [121 spectra from bin A 
and B, the regions of highest excess flux, and spectra from bin C 
and D, where little excess flux is seen. 

3.5 Velocity of the hot gas 

Assuming a Gaussian profile for the [Fe X] A6374 emission line de- 
tected in box 1 we can fit the velocity width of the detected broad 
excess flux. We do this for each of the methods of continuum sub- 



traction, as described above. The FWHM velocity width of the fea- 
ture determined from the different continuum subtraction methods 
is 300 to 650 kins" 1 , from box 1 (the velocity width of the [O 
I] A6363 is ~ 200 km s _1 ). It should be noted that the predicted in- 
tegrated line profile of the coronal line emission is model dependent 
dSarazin & Granevl 1 99 ll) , however our data are not deep enough to 
allow for an investigation of the line profile. The thermal line width 
of the coronal line emission is typically only 20 — 30 kins" 1 , so 
the broad w idth of this feature is like ly due to turbulent motions of 
the hot gas. lSarazin & Granevl dl99lh show that for a gas of 2 mil- 
lion K with maximal turbulent broadening (velocities limited by the 
sound speed in the hot ambient medium) predicted velocity widths 
at FWHM are approximately 1700 km s _1 . The motions in the gas 
we observe are therefore highly subsonic. 
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Box 1.° 


Th 1 1 1 v ( over fm 

Box 2. a 


2 -1^ 

fe ) 


All pixels. b 


[FeXIV]A5303 


< 1.9 x icr 16 


< 1.7 x 10" 


-16 


< 4.9 x 10- 16 


[CaXV]A5445 


< 4.4 x 10- 17 


< 3.7 x 10" 


-17 


< 3.9 x 10- 16 


[CaXV]A5694 


< 1.6 x 10- 16 


< 5.9 x 10" 


-17 


< 1.4 x 10- 16 


[FeX]A6374 


1.5 x icr 16 ± 2.4 x icr 17 


< 3.0 x 10" 


-16 


1.1 x 10~ 15 ±3.5 x 10~ 16 


[Ni XV]A6702 


< 4.6 x 10- 17 


< 6.6 x 10 


-18 


< 6.0 x 10- 17 



Table 2. Same as Column 1-4 of table [T]in units of flux. 



4 DISCUSSION 

We have conducted a deep search for 5 species of optical coro- 
nal line emission, specifically [Fe XIV]A5303, [Ca XV]A5445, [Ca 
XV] A5694, [Fe X] A6374 and [Ni XV] A6702. These lines probe gas 
at temperatures of 2, 5, 5, 1.3 and 2.5 million K respectively. 

We report four upper limits and a 6.3 sigma detection of [Fe 
X]A6374 emission. This detected feature is broader than the Ha 
nebulosity at 10 4 K, having a FWHM velocity width of about 
300—650 km s -1 and is blueshifted with respect to this cooler 
emission. [Fe X]A6374 is emitted from gas at one million K, the 
lowest temperature of the coronal lines in our wavelength range. 

Assuming that coronal line emission is a product of the low 
level of residual cooling out of the hot gas observed in the central 
galaxies of galaxy clusters we can determine the rate of cooling 
from the flux of these lines. Upper limits for the detection of [Fe 
XIV]A5303 and [Ni XV]A6702 are consistent with t he X-ray cool- 
ing ra te, in the absence of heating, of — 10 M yr -1 JSanders et al.1 
2008) above 0.8 keV, over a region 20 arcseconds in diameter. 
However in the case of [Fe XIV]A5303 we stress the difficulties in- 
volved in the subtraction of the stellar continuum which introduces 
a larger error into the detection limit (see section [37TI and Fig.[9]). 

The rate of cooling implied by the detection of [Fe X]A6374 
emission is a rate of — 20 M© yr -1 in the central 20 arcsec 2 , twice 
that inferred from X-ray observations and from the upper limits 
of the other coronal lines (see table Q]). If this emission is solely 
due to the cooling of the ICM in this cluster we would expect 
similar or greater detections from the [Fe XIV]A5303 and [Ca XV] 
lines. An alternative for the iron emission is that it originates in gas 
which was cooler and has been heated. There are lots of possible 
sources of heating including shocks, photoionisation from central 
active galactic nuclei (AGN), thermal conduction, turbulent mixing 
or heating associated with the radio bubbles. We caution that re- 
cent changes in the atomic parameters indicate that the emissivity 
of [Fe X]A6374 is larger than previously thought. This would im- 
ply a slightly lower cooling rate than calculated here and will be 
investigated in further work by Ferland et al. (in prep). 

C IV 1549 A emission has been discovered by ISparks et"aD 
(120091) in the outer regions of M87, co-located with its Ha fila- 
ments. This probes the slightly lower temperatures of — 10 5 K. O 
VI 1032, 1035 A emission (— 1 55 K) has also been seen in a small 
samp l e of cool core clusters (lOegerle et al.l l200ll ; iBregman et al.l 
l2006h . iBregmanet all d20Q6h note that in Abell 1795 the cooling 
rate implied by the 10 5 5 K is larger than that implied by the X-ray 
emission, as is the case here. They suggests that non-steady cooling 
of material may be the cause of this conflict. 

Most AGN show emission from forbidden, high-ionisation 
lines, the so called extended coronal line regions (CLRs) 
(iMullanev et al.ll2009l : IMazzalav et al.ll2010h . This emission is be- 
lieved to be from a region outside of the broad-line region and 
inside the narrow-line region and lines are often observed to be 



blueshifted by 500 or more km s" -1 . IMazzalav et all bOlCh find 
that photoionisation from the central AGN is the likely major ion- 
isation mechanism for these regions. Their results also confirm the 
observation that higher-ionisation lines are emitted from a more 
compact area while the lower-ionisation lines can be found further 
into the narrow-line region. 

For a typical Sey fert galaxy with an ionising luminosity 
Li on = 10 43 5 erg s -1 JFerguson etaD dl997h find that the region 
containing [Fe X]A6374 emission should be restricted to the inner 
20 pc. The size of the CLR scales with the ionising luminosity of 
the central source as L\{2 JFerguson et al.|[l997l) . Box 1, where we 
find significant [Fe X]A6374 emission, is located — 15" from the 
centre of the r adio source, a distanc e of —3 kpc. Scaling up from 
the values of JFerguson et al.l lT997), a region of the size —3 kpc 
implie s a huge ionising luminosity of 10 47 ergs -1 . iTavlor et al.1 
(120060 have found the nucleus of NGC 4696 to be very faint. Their 
upper limit on the luminosity from X-ray observations between 
0.1 — 10 keV is — 10 40 ergs -1 . The spatial extent over which 
we find the coronal line emission leads us to rule out the possibility 
of photoionisation by the central AGN. It is possible, looking at the 
residual flux over the spectral region where we would expect [Fe 
X]A6374 emission, that the spatial distribution of the 1 million K 
gas (Fig.[TTJ traces that of the 10 4 K filaments and large dust lane 
however this is a tentative result and would need confirmation with 
deeper observations. 

Another explanation for the excess flux such as emission from 
a merging clump of gas or another filament is unlikely, as any ex- 
tra component of [O I]A6363 emission should also be present in 
the [O I]A6300 emission and no red wing is observed on this line. 
One explanation for the blueshift could be that the radio source is 
interacting with and driving the plasma away possibly entraining 
with it cooler gas. Fig. [T3l shows a two component fit to the [N II] 
and Ha emission in the box 1 spectrum. There is evidence for a 
blueshifted component to these lines and the best fit gaussians im- 
ply a blueshifted velocity of — 140 km s -1 . This is less than the 
observed blueshift of the [Fe X]A6374 line (-230 km s -1 ) so is 
inconclusive. 

Sanders ctalJj2008h have used deep XMM-Newton Reflection 
Grating Spectrometer (RGS) observations to show that the centre 
of the Centaurus cluster contains X-ray emitting gas down to at 
least 4 million K through detections of the Fe XVII line and lim- 
its on the O VII emission. The metal abundances are ve ry high in 
the inner 30 kpc of this cluster JSanders & Fabianl2006al) . The iron 
abundance of the hot gas is —2 times solar and the calcium abun- 
dance —3 times solar. The Chandra spectra show evidence of an 
off-centre abundance peak, at —20 kpc. The best traced elements 
of Fe and Si then show a decline in abundance in the nucleus. 

At these high abundance values we expect to detect, in our op- 
tical spectra, lines of [Ca XV], probing gas of 5 million K, detected 
already by its X-ray Fe XVII emission. Contrary to this expectation 
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Figure 11. The spatial distribution of the 1 million K gas. A bove: The [N 
IIJA6 583 emission, overlaid with VLA 5 GHz radio contours ( iTavlor et"al] 
120071 . The white box indicates the region where we have our deepest ob- 
servations. Below: This region binned on a 5 by 5 fibre basis showing the 
residual flux at the position of the [Fe XJA6374 emission across the field of 
view. The brightest bins are regions of highest flux. The spectra from bins 
A-D are shown in Fig.^1 For comparison box 1 and 2 are marked in red. 



we do not detect any significant calcium emission. Upper limits in 
most cases suggest the rate of radiative cooling is less than that 
inferred from the X-ray spectra of hotter gas and provide a more 
stringent constraint than the other coronal emission lines. 

We explore four possible explanations for this result; inaccu- 
rate continuum subtraction; a lower cooling rate; a low calcium 
abundance and calcium deficiency in the intermediate temperature 
gas. We will deal with each of these points in turn. 

The continuum subtraction has been attempted in a variety of 
ways (see section 13.11 ) and in all cases produces similar results. 
Where these results were upper limits and not clear detections the 
largest value was taken and is given in table Q] Where the stellar 
spectra exhibit many features the uncertainties are large and the 
detection of broad line emission is very difficult. However, in the 
absence of a template spectrum in which we are sure there is no 
emission from intermediate temperature gas and observations with 
similar abundance ratios across the field of view this is the best we 
can do. 

The cooling rates are derived for cooling without heating and 
so may be lower than stated, they may also vary across the tern- 



Figure 12. A comparison of the spectra from four of the regions marked 
on Fig.[TT] The spectra have been continuum subtracted and the [O I] A6363 
object emission line and the two sky lines to the left of the Oxygen line have 
been fitted and removed. The bottom panel shows the residuals between the 
two spectra. 



perature range probed. Investigations of the quantity of molecular 
gas and dust and of recent star formation in these objects may help 
constrain the quantity of cooling, to very low temperatures, that 
we expect. A much lower cooling rate may provide an explana- 
tion for the lack of [Fe XIV] A5303 and [Ni XV] A6702 emission we 
see however the calcium lines should originate from a temperature 
range that has been directly observed from measurements of the Fe 
XVII line ratios, so the lack of this emission is not easily under- 
stood. 

The Centaurus cluster has very high metallicities which 
peak at a radius of 20 kpc, then decline towards the nucleus 
JSanders & Fabianll2006ah . This enhancement and decline in the 
metallicities are reproduced using both XMM-Newton, Chandra 
CCD observations and XMM-Newton RGS observations. An off- 
centre abundance peak with a depression in abundance towards 
the centre of clusters and groups is often observed in objects 
with cool cores (see for e xa mple | De Grandi & Molendi _ 200 ll ; 



iRasmussen & Ponmanll2007l and lDe Grandi & Molendill2009h 

There are a number of biases and uncertainties in abundance 
measurements from spectral modelling. Inaccuracies in modelling 
the te mperature st ructure of galaxy clusters can introduce an 'Fe 
bias' (Buote 2 000h . Here the measured Fe abundance is lower than 
the actual value due to fitting a multiphase gas with single or only a 
couple of temperature components. The opposite effect, an 'inverse 
Fe bias', is also seen in cool core clusters when a single temperature 
model is fitt ed to a multiphase plasma with temperatures between 
2-4 keV feasia et al.l l2008). In the case of Centaurus the high 
spectral resolution of XMM-Newton RGS allows the temperature to 
be constrained with line strength ratios giving a more robust chec k 
on the temperature components of the models dSanders et al.l2008b . 

Sedimentatio n in the centre of galaxy clusters 

jFabian & Pringld Il977l) can cause the metal abundances to 
rise, and could be reversed by the effects of thermal diffusion 
dShtvkovskiv & Gilfanovl 120101) . Neglecting the effects of reso- 
nance scattering also underestimates the abundances of metals 
in clusters. This effect is at most 10 per cent, so cannot fully 

explain the central abundance dips observed in galaxy clusters 
1 — 1 ■ 11 ■ — 1 

(Sanders & Fabian 2006b). The drop in abundances in the Centau- 
rus cluster may be explained with a complex model involving three 
temperature components and additional absorption however the 
errors on the inner-most radial bins become s ufficiently large tha t 
it is impossible to tell if the drop is real or not dFabian etaU2005l) . 
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Figure 13. A two component fit to the Ha; and [N II] emission from box 1 . 
There is a evidence for a slightly blueshifted component to these lines with 
a blueshifted velocity of ~140kms — 1 . This is smaller than the velocity 
shift seen in the [Fe X]A6374 emission line of ~230 km s _1 . 



A major contributor to a central abundance drop must be de- 
pletion on dust. The central few kpc of the hot gas will be domi- 
nated by stellar mass loss, in which most metals are bound in dust 
grains. They will slowly be introduced i nto the hot phase by sput- 
tering, at a ra te dependent on grain size dPraine & Salpeterj|l979l : 
lBarlowl[l978h . Much iron however could be injected into the hot 
phase through SNIa. 

Calcium is one of the most depleted of all refractory ele- 
ments in the presence of dust <lFieldlll974l : rSpitzer & Jenkin s Hl975l : 
Cowie & Sonsaila with depletions of 10 4 relative to so- 

lar abundance typical in dense clouds. NGC 4696 hosts a huge 
dust lane which almost completely encircles the core and spirals 
out to the north east, tracing the morphology of the Ha filaments. 
The X- ray absorption colum n density is also highest in the same 
region dCrawford et al.ll2005h . It has moreover a large quantity of 
infrared-emitting dust in its nucleus (jGoudfrooii & de Jon3 1 19951 ; 
iKaneda et al.l2007l) . lKaneda et all ( 120071) show from Spitzer Multi- 
band Imaging Photometer (MIPS) 24 /am surface brightness pro- 
files that the dust emission increases steeply in the inner 30 arcsec- 
onds (6.3 kpc), a region slightly larger than our field of view and 
where we observe the metallicity drop. 

The cumulative gas mass obtained f rom th e X-ray emitting hy- 
drogen density profile of iGraham et al.l d2Q06h in the inner 30 arc- 
seconds or 6.3 kpc (10 arcseconds, 2.1 kpc) is 9.4 x 10 8 M 
9.2xl0 7 M ) Assuming a dust to gas ratio similar to our Galaxy 
Crawford et al1l2005l : lEdge et alJl201Qh we get a total dust mass 
of 9.4 x 10 6 M (9.2 x 10^ M ). The dust mass estimated 
bylGoudfrooij & de Jond dl995h using the Infrared Astronomical 
Satellite (IRAS) is 4 .7 x 10 6 M and the Spitzer MIPS result from 
IKaneda et all (120071) is 4.2 x 10 5 M . These values are considered 
to be lower limits to the dust mass as the instruments are insensi- 
tive to very cold dust and imply that the inner interstellar medium 
in NGC 4696 is highly deficient in refractory elements. 

A number of cool core clusters have been fround to be defi- 
cient in calcium in the warm (10 4 K) emission line nebulae, where 
[Ca II] and [Ca V] emis sion lines would be expected (iFerlandl 19931 : 
lDonahue& Voitlll993l) . Our spectra of NGC 4696 are consistent 
with this since they show no evidence for [Ca II] or [Ca V] emis- 
sion in the inner regions of NGC 4696. 

These results indicate that NGC 4696 is deficient in gas-phase 
calcium at temperatures below ^5 million K. The lack of this warm 



and intermediate temperature gas phase in the inner 20 arcseconds 
implies the calcium in dust has never been part of the hot ICM and 
was probably introduced in dusty stellar mass-loss. 

The high calcium ab undance in the hot gas, discovered by 
I Sanders & Fabianl (l2006al) , has its greatest contribution from the 
hot 4 keV emission (the contribution from the best fit to the lower 
temperature component, 0.5 keV, is negligible). This is consis- 
tent with the dust being sputtered by the hotter gas at large radii 
(^few kpc), where there are also no shielding 'cold' filaments. 



5 CONCLUSION 

We report the detection of [Fe X]A6374 coronal line emission in 
NGC 4696. This emission probes temperatures of 1 million K and 
is the first detection of intermediate temperature gas in this object. 
We fail to detect emission from coronal lines of higher temperature 
gas including those of [Ca XV] which probe gas of 5 million K 
and which we had expected to detect due to the high abundance of 
calcium in the hot X-ray emitting gas. 

We conclude that calcium is likely to be depleted in the dusty 
central regions of NGC 4696. This is consistent with our apparent 
lack of [Ca V] and [Ca II] ions which probe lower temperature gas 
and with the negligible contribution to the abundance in the hot 
gas of the lower temperature 0.5 keV calcium lines. The dust in the 
central region of the galaxy is likely due to stellar mass-loss and has 
survived as dust grains due to the shelter of surrounding cooler gas. 
The abundance of calcium is higher in the outer 4 keV ICM since 
dust is sputtered there where the dust is less protected by the cooler 
surroundings and the cold filaments. Deeper, high resolution X-ray 
observations which better constrain the central abundances of Fe, 
Si, Ca and Ne would help to distinguish between the processes that 
contribute to the metal abundance. 

There is now strong evidence to show that the central galaxies 
in many cool core clusters are pl aying host to large quantities of 
dust and cool ga s (see for example [McNamara & O Connellll 19921 ; 
lEdge et al.l2010l) . lDonahue & Voitl dl993h have also shown that the 
[Ca II] doublet is much weaker than expected in a sample of BCGs. 
They conclude similarly that calcium is most likely depleted into 
dust grains in their sample. 

The cooling rate inferred by the [Fe X] A6374 emission is large 
(20 M yr _1 in a spatial region of 20 arcseconds 2 ) compared with 
that determined from other lines in the optical and X-ray spectrum. 
This and the apparent lack of [Fe XIV]A5303 and [Ni XV]A6702 
emission rule out a steady cooling flow from the 10 7 K gas in this 
object. Some gas may however be cooling non-radiatively for ex- 
ample by mixing with the colder gas. The strength of the [Fe X] 
emission suggests that the million K gas is being heated rather than 
condensing out of the hot ISM. 
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